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The origin of the lightest scalar mesons is studied in the framework of instanton liquid model 
(ILM) of the QCD vacuum. The impact of vacuum excitations on the a-meson features is 
Q I analyzed in detail. In particular, it is noticed that the changes produced in the scalar sector 

I may unexpectedly become quite considerable in spite of insignificant values of corrections to the 

ly^ ' dynamical quark masses and then the medley of cr-meson and those excitations may reveal itself 

as broad resonance states of vitally different masses. 



Nowadays well-known theoretical results of existing scalar mesons, the non-zero vacuum expec- 



> 

m 

■ tation value of which is strongly argued by chiral symmetry breaking, find earnest experimental 
^ support. It comes from studying the low energy S'-wave of vr-vr-scattering where the presence of Solt- 
ys tary low mass scalar resonance looks inevitable pQ and the radiative 0-meson and heavy quarkonia 
O I decays as seen from the view point of chiral shielding idea [2]. Thus, the present situation in the 
^ • subject entirely justifies the theoretical expectation that the physics of scalar mesons is driven by 

the Goldstone dynamics [3J. 
P-i! However, it is still difficult to understand the properties of scalar mesons (resonances) in terms 
r-| of the QCD basic fields, in particular, the quark and gluon origin of light a-meson (/o(400 - 1000 
I MeV)) and broad resonance /o(1000 - 1600 MeV) [i]. In a sense, the current status of the lightest 
scalars could be summarized by the strange assertion that people know where the scalars are but 
do not know what they are. Meanwhile, the solution of this still pending problem can be found 
by studying the QCD vacuum structure [5] and searching putative quark-gluon plasma state and in 
recent time the critical point in ultra- relativistic heavy ion collisions [B] . An idea to exploit cr-meson 
(precisely, its coupling to photons) as a key guide to exploring the existence and features of mixed 
phase of strongly interacting matter created in nuclei-nuclei collisions at lower energies [7] requires 
such an insight as well. Attempts to find an explanation of splitting two lowest scalar mesons in the 
framework of rather sophisticated models mixing the quark-antiquark states with the glueballs (see, 
for example [8],[9]) are not fully successful and the results are ambiguous. 

In the present paper we consider the origin of lowest scalar mesons in the instanton liquid model 
(ILM) of the QCD vacuum [TU] mixing its phonon-like excitations with the scalar mesons treated 
in a standard way as bound quark-antiquark states. The effective Lagrangian of phonon-like excita- 
tions is similar to the dilaton one [8] and describes the state with the quantum numbers J^*" = O"*"^ 
(J = total angular momentum, P = parity, C = charge-conjugation eigenvalue) of glueball. One 
specific feature of this Lagrangian is the form of its kinetic term and rather strong interactions with 
quarks (resulting in the strong mixing with the cr-meson field p2j). 

Avoiding the technical details of constructing this effective Lagrangian we mention briefly here the 
major steps necessary for the developed approach. In particular, it is grounded on the hypothesis the 
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vacuum field configurations are stabilized at the certain characteristic scale and their action develops 
a well-defined minimum at the point of average configuration size. The generating functional in ILM 



D[A] e-^(-^) , 

where >S'(^) is the Yang-Mills action, is supposed to be saturated by instanton superposition 

N 

-A;:(a:) = E^;:(^;7.). (i) 



Here A^(a;;7j) is the (anti-) instanton field in the singular gauge 

A» = ^ u;^,,. ^ , y = x-z, (2) 

where ji = {pi,Zi,uJi) are the parameters characterizing the i-th (anti-) instanton of size p with a 
matrix of colour orientation u and with coordinates of its center position at z and g is the strong 
coupling constant. For anti-instanton the 't Hooft tensor should be substituted according to 77 — >■ r/. 
We do not discuss here the mechanism of instanton ensemble stabilization at the scale of p (see, for 
example p3]) but focus on extracting some phenomenological results in the context of our interest. 

The form of action with peculiar minimum for the saturating configuration equilibrated makes an 
existence of "oscillations" (one should keep in mind we are working in the Euclidean space) of this 
configuration around the very natural size p. In principle, such a description should be done by the 
corresponding Green function. However, calculating it is not simple and eventually impracticable 
because the Green function is found to be singular [14j. It was noticed in Ref.[15| that if one is inter- 
ested in the major terms of a generating functional exponent then less information on interrelation 

between the saturating configuration A and field of "oscillations" ^ at the characteristic scale p is 
necessary to calculate the respective kinetic term of the effective Lagrangian. The instanton ensemble 
excitations generated by a certain impact might naturally be described by more general saturating 
configurations of the form of Eq.([2]) if the instanton size and its colour space orientation are varied 
(p — R{x, z), uj""^ —>■ Q^-^'lx, z)) and the deformation fields could be defined by dealing with minimal 
action requirement 6S = 0. As a result it leads to a more accurate (than plain superposition Eq.([l])) 
solution and allows to determine the interrelation looked for. In fact, it can be carried out directly 
since the deformations are defined by the multipole expansion done in the center point of instanton 

Rin{x,z) = p + c^y^, + c^^ y^,yu + ■■■ , \y\ < L 

y y Vi^ 

Rout{x,z) = p + + d^u ^ + ■ ■ ■ , \y\>L, 

y y y 

(similar expansion should be done for instanton orientation in the "isotopic" space Q{x,z)), here L 
is a parameter which fixes a sphere radius where the increasing multipole expansion with distance 
increase changes its behaviour for the decreasing one because of the requirement of deformation 

regularity. The coefficient in Eq.® just corresponds to the function which we are interested 

in (by the way, for the sample of deformation field Eq.Q it is valid ~ ~Sz~ ^^"^ interrelation 

between two fields is defined by the solution of Eq.([2]) with R{x, z) included). 

Now calculating an action of this crumpled (as was called in Ref . pjj| ) configuration we are able 
to determine the kinetic term of effective Lagrangian at the characteristic scale p (certainly, we can 
not profess a higher precision in this quasi-classical approximation as it was mentioned before) in 
the form 

S,,^ = j dx\G1M)GlM)-p = ^{5,p)\ ^=^f^. (4) 
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here (3 = Sn'^/g'^ is the single (anti-)instanton action at the p scale which in its general form can be 
presented as 

sip) = Pip) + 5 In(Ap) - In ^^^^ + np^p^ , (5) 
with the function Pip) = —IhCn^ — 61n(Ap) , A = Ajjg = 0.92Apy. , and the constant Cjy^ 
depending on the renormalization scheme Ctv^ ~ ^n^N^i)\{N^-^2y. ' ^ ~ ^ ^ b = — — ^-g— — 
Nf is the number of flavours, Nc is the number of colours and P = Pip), P = P + IuCa?^ are the 
magnitudes of Pip) function at the fixed value of p, ^ is a constant characterizing the repulsive 
power of pseudoparticles = xt^T^^ ^ ^'^^ instanton liquid density. Holding the terms of 



second order in the small deviations from the point of action minimum 



dsip) 
dp 



we receive the 



approximate result 



sip) ^ sip) + ^^ (6) 



where s*^^)(p) ~ ^ ^^^^ 

deviations from the equillibrium value Pc = P y- ~ — P- Finally the deformation field could 

be described by the following effective Lagrangian density [TT].|15j 



=, and the scalar field ^ = bp = p — pc — p — p just realizes the field of 
p2 

1/2 



TiK I / 
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^. = ^<i V-g) > (7) 

with the mass gap of phonon-like excitations as = = 4= (the coefficient k = 0.9/3 

generates the scale of order 1 GeV with A ~ 280 MeV). 

Then the quark determinant Z^p for the stochastic ensemble of pseudo-particles looks like 



where Si^p, A) is the action of QCD with massless quarks which should describe the spontaneous 
breaking of chiral symmetry and generation of dynamical quark mass of order 300 MeV [I6] . Eventu- 
ally the quark generating functional (after absorbing the variations of (anti-)instanton average size) 
takes the following form (for Nf = 2) of the integral over saddle point parameters [12] 



Z^ = J dX DM D-^^Dip e ^^^^ , 

N r ^— 

= -N\n— + 2 j dx iNf-l) \ ""f-' iDetM) 



N 

/V In- 

^}(/c) (27r)'^5(A; - /) (^-k + imfg vik, A;)) +imfg uik, I) (pik - I) ^g(/) , 

where is the total number of pseudo-particles in the volume V in = N/V), M is the Nf x Nf- 
matrix of scalar boson fields in the fiavour space and /, g are the fiavour indecies. The vertex functions 
vik, k) and uik, I) are defined by the zero-modes (those are the solutions of the Dirac equation in 
the field of (anti-) instanton of zero energy) and have the following forms 

vik, k) = G\k), uik, I) = = Gik)G'il) + G'ik)Gil) , 

dp 



moreover 



Gik,p) = 27rpF(fcp/2), G'ik,p) = ^Gi^A^ ^here 



dp 

Fix) = 2x [Ioix)K,ix) - hix)Koix)] - 2 hix)K,ix) 
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and Jj, Ki {i = 0, 1) are the modified Bessel functions. 

After executing tlie bosonisation procedure tlie Lagrangian density takes tlie following form 



^ — ^ip~^ ^■4>,ip ' (9) 

which describes the quarks and phonon-like excitations of instanton liquid. The Lagrangian param- 
eters are fixed by saddle point available as 

d£_ a/: _ 

9A~ ' dM~ 

In the particular case of negligible impact of field saddle point values of parameters are fixed by 
the requirements 

1 ^— N 

"J (27r)4 k'' + m)v\k) V ' ^ ' 

here v{k) = v{k, k). 

As it was mentioned above the Lagrangian density ([9]) is quite similar to the dilaton one. How- 
ever, there are, at least, two substantial distinctions. First, it is the presence of factor uk in 
the kinetic term. This factor is analogous to the factor of chiral vr-meson Lagrangian density 

£^ = -ij- dfj,Trd^Ti, and sets up the scale of phonon-like excitations. Average instanton size p = pf) 
is quite relevant dimensional unit for phonon-like field and, hence, the "correct" dimension for the 
factor in kinetic term could be 

= nnf . 

The second distinction concerns the interrelation of phonon-like field and quarks which does not 
disappear in the chiral limit and is defined by the vertex function u{k, I). 

Then a nontrivial solution of saddle point equation Eq.([Tp) fixes the dynamical quark mass 
M(k) = mfv{k) together with the quark condensate 

•/ ,t M ..r f dk Mik) 

The excitations of quark condensate with the respective quantum numbers are seen as a- and vr- 
mesons (in this paper we are dealing with Nf = 2 only) and 

M = m (1 + a) e*^""" , 

where are the Pauli matrices. The corresponding correlation functions -Ro- and Rj^ are given, as 
known, in the form 



f dki Mf f dki ((kik2) + M1M2) M1M2 , , 

( \ -4 OAT AT f {{klk2) - M,M2) M^M^ 

with k2 = ki + p, Ml = M{ki), M2 = M{k2)- Expanding the integrals around small values of 
momentum p as 

Rn{p) = Pn P'^ , Ra{p) = + Pa p"^ , 
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we are able to calculate the following decay constants 

r dk M^ik)-^M'ik)M{k) + ^{M'{k)r 
2 J (27r)4 {e + M\k)y 

^ r dk \ 2M\k) ^ M\k) 



{2ny 



Ai- 



{k^ + M\k)y ' (k^ + M^k)) 



Ai = — M(k)M'(k) - M'(k)M'(k) + M(k)M"(k)] , 
16 \k / 

^ _ k^ + M^{k) ^ M{k)M'{k) ^^ k^ ^ M{k)M'{k) \j ^ 
k^ + M^k) ,./,^,.///,^ M(k)M'(k)\ 

H g ^ ^ iM'{k)M\k) + M{k)M"{k) ~k^] ' 

(here the prime denotes the derivative calculation) and the cr-meson mass as well 

dk - M^{k)) M^{k) 



(27r)4 {k'' + M\k)f ■ 

The width r^. of cr-meson decay into tttt is defined by the set of graphs a — 5 (see Fig.l) for the 
anharmonic term ~ cttt^ of effective Lagrangian 

f dmdl r dk ,„ ^, .. ^n,,. ~n, 

'""^ ^ J (2^ J (2^ W^{k,m,l) a{m) 7r"(/) 7r"(m - /) , 

W^{k,m,l) M^{k) 



AN^ ~ A;2 + M\k) 

k{k - m) - M{k)M{k - m) 
{k^ + M\k)){{k - mf + W^ik - m)) 

{k^ + M^{k)){{k-lf + M^{k-l)) ^ ^ ^ ^ 



M{k)M{k - m) + 



k{k - m)M{k -I)- k{k - l)M{k ~ m) - {{k - m){k - l)M{k) - M{k)M{k - m)M{k - I) 
^ (fc2 + M\k)){{k - mf + W^ik - m)){{k - If + M^{k - I)) ^ 

xM{k)M{k - m)M{k - I) . 

(four terms of this expression are in explicit correspondence with the four graphs of Fig.l). 




Figure 1: The single line corresponds to the quark contribution, the double one corresponds to 
(T-meson, the short dashed line shows the phonon-like field and the long dashed one shows 7r-meson. 

Then for the decay width of cr-meson in its rest frame we find that after integrating the function 
Wa{k, m, l)/Fcr/F^ over the loop at fixed momenta of vr-meson [17] we have 

8nM^ Ml) ' 

here denotes the integral over anharmonic contribution of type aii^ . At zero outgoing 
momenta the anharmonic contribution disappears VFcr(^;0,0) = because in the chiral limit vr- 
meson fields may appear in the combination with the derivative term (9^vr)^ only. As known, this 
fact, in some extent, results in the small cr-meson width. 

In the meantime, if the phonon-like excitations are not taken into account tt and cr-mesons are 
defined by the following fundamental characteristics 

~ 121 MeV , ^ 96 MeV , 
~ 550 MeV , V„ ^ 2.8 MeV , 

and TT-meson mass can be extracted from Gell-Mann-Oakes-Renner (GMOR) relation. The saddle- 
point parameter m taken as dimensionless quantity is equal to m = 4.817- 10~^ . Then we have for the 
instanton liquid density nj = 1.03, for the average pseudo-particle size pA = 0.28 and eventually 
for single instanton action at this characteristic scale /3(p) = 19. The dynamical quark mass and 
quark condensate read as 

M ~ 385 MeV , - ^ -(381)^ MeV^ , 

at A = 280 MeV. Surprisingly, the small width of a-meson is a result of interference of all terms 
contributing. In particular, in dimensionless units the corresponding contributions are as follows 

a = -1.199-10-2, /5 = -7.137- 10-^ 7 = 2.066-10-2, 5 = -1.695- 10-^ a+/5+7+5 = -1.648- 10-^ . 

Actually, many quantitative estimates are based on the calculations of 6 graph only and it results 
in the magnitude of width T^r which is declared to be a few hundred MeV. Assuredly, there is an 
unpleasant drawback of this calculation. It comes from the opposite signs of contributions of a and 
7 diagrams because the obvious compensation produces rather small value which is poorly controlled 
by one-loop calculations. Unfortunately, for the time being two-loop calculations are hardly realized 
technically. 

In order to include the changes coming from phonon-like fields entering the game, the perturbation 
scheme was proposed in Ref.[T2] in the tadpole approximation. Actually, it becomes possible, because 
of the quark condensate presence in the following sum 

to keep the first term only in the respective inserts into quark Green functions and vertices. Then the 
dynamical quark mass is defined not only by the vertex v but the tadpole contribution with vertex 
u as well (see Eq.(j8])). Amazingly, the dependence of diagram contributions on the coefficient k is 
cancelled in this approximation. Generally, there is hardly any impact on the instanton liquid (it 
is negligible), average instanton size and parameter /3(p) do not change but instanton liquid density 
slightly increases n/A^ = 1.11. The saddle point parameter is getting smaller m = 3.481-10-^ which 
results in the dynamical quark mass and quark condensate decreasing (see Fig. 2). 

M ~ 324 MeV , - ^ -(343)^ MeV^ . 
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Figure 2: The dynamical quark mass as a function of momentum. The dashed hne shows that 
behaviour when the phonon-hke fields are not taken into account. 



However this quite noticeable decrease of quark dynamical mass (and quark condensate as well) 
should not be taken very seriously because the scale of A in the ILM is not strictly fixed and it 
is normalized by calculating the observables. Thus, if one considers, for example, vr-meson there is 
only the possibility to fit parameters (for both with phonon-like excitations of instanton liquid and 
without them) in such a way to get the satisfactory value of constant F^, the GMOR relation and 
etc. Then characteristic parameters of field scales in the tadpole approximation are the following 

~ 106 MeV , ~ 86 MeV , ~ 337 MeV . 

Peculiarly, the field (f develops remarkably larger scale comparing to tt- and cr-mesons. 

The width of phonon-like field (/? decay into vrvr-mode is defined by the diagrams similar to 
a and 7 but with changing cr-mcson line for the line of phonon-like field. Then the corresponding 
anharmonic term of effective Lagrangian density <^7f^ takes the form as 

" y (2^y (2^ W-^(A;,m,0 <p{m) 7:\l) 7r"(m - I) , 



W^{k,m,l) _ M{k) AM{k) 
iTVc ~ k^ + M^{k) ^ 

k(k-l)+M{k) M{k-l) ,,,, , nA,./, , ,N 

two terms of the sum here correspond to the diagrams mentioned and AM{k) = ruf u{k, k). Finally, 
we have 

3 / 4M2\'/' 



where W^p looks similar to but with changing the factor in denominator providing the dimension- 
less result for F^F^). The mass of phonon-like field and its width are 

~ 827 MeV , ~ 340 MeV , 
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but these characteristics for a-meson become as 

~ 1.004 MeV , ~ 270 MeV , 

(for quark loop without the tadpole contribution the width of phonon-like excitation is estimated to 
be r<^ ~ 780 MeV). The partial contributions of diagrams to the a-meson width read as 



a 



-7.508-10^3^ (3 = -2.205-10-3, 7 = 8.567-10-^, 6 = -2.999-10-\ a+(3+-f+6 = -1.447-10-3. 



It is instructive to compare the magnitudes of a^j and jS^j integrals contributing dominantly to the 
(T-meson mass (upper line corresponds to the contribution without phonon-like field) 



a,, = 1.333-10-3, = 1.812- 10-^ = 2.872- 10-^ , 

= 3.508-10-3, = 1.427-10-2, = 2.202 -lO'^ , 

(for all that the packing fraction parameters are np^ = 5.997-10-3 for the upper line and np^ = 
6.494-10-3 for the lower one). As seen, the coefficient a^^ gets the strongest impact. 

Now let us turn to the problem of mixing two scalar fields if their Lagrangian density tolerating 
the decay into two 7r-mesons mode has the following form 

-£ = -S. + M^) ^'-^ + AO - 

__iL (fc2 + M^) n'^ + A^a + W^ ^n"" + W„ _ 

Here the parameter regulating the component mixing A is given by the tadpole graph of Fig.l (and 
the sign of mixing term coincides with the sign of quark condensate) and could be presented as 

As it was mentioned above, unlike the dilaton Lagrangian [Hj in this scheme the field mixing survives in 
the chiral limit and for the instanton liquid parameters used above we have A = = (483 MeV) 2. 

The masses of diagonal fields are given by 

, Ml^Ml ((M2 + M2)2 + 4A2)V2 

Then the minimal mass of light scalar component is limited by constraint M_ = 2Mn-, otherwise it 
could be stable in strongly interacting mode, and the upper limit of the A magnitude responsible for 
mixing looks like 

what for existing parameters gives A^ax = (868 MeV) 2. Comparing to the standard A value one 
may conclude the mixing effect is quite significant. 

Passing through the standard procedure with 6* as a mixing angle 

(f = cos 9 ip' + sin 9 a' , 
cr = — sin 6^ + cos 9 a' , 

we receive for corresponding widths of and o' the following results 

W^. = W^cos9 - W^sm9 , W^^ = W^^sin^ + W^cos9 , 
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where the value of mixing angle could be obtained from 



It becomes clear from Eq. fll3l) the mass of lighter components is getting smaller whereas the heavier 
one gains more mass. Remembering the phonon-like excitation was lighter before mixing we find out 
for if' the following characteristics 

M^/ ~ 749 MeV , T^, ~ 265 MeV . 

The heavier component a' becomes as 

~ 1.063 MeV , ^ 175 MeV . 

The mixing angle is ^ = 27.5 degrees. Unfortunately, the precision of calculating the particle decay 
widths is not high and allows to rely on the order of magnitude obtained only. The reason is 
tightly linked to the contributions of the a and 7 graphs which are considerably larger than others. 
Calculating the a and 7 contributions within 10 % precision gains a much poorer precision for the 
total result. In general it is rooted in that fact the estimates are obtained at the limit (sometimes 
beyond it in other approaches) of applicability of expansions used. 

Summarizing, we would like to underline the approach proposed here sheds more light on the 
possible origin of scalar mesons. Assuredly, many other models have expressed the similar intends 
and succeeded to different extents. However, the calculations carried out in this paper demonstrate 
a quite robust effect of phonon-like excitations of instanton liquid in the scalar sector and capacity of 
such a model to successfully absorb new results and hints of running experiments. It concerns, first 
of all, the meson state splitting and its strong coupling scale. For example, both quark-antiquark 
and glueball O"*""*" operators mixed to create the scalar state in lattice QCD calculations give the 
mass of the lightest meson essentially suppressed with respect to the mass of corresponding (0^+) 
glueball and demonstrate strong lattice spacing dependence |18|. Besides, the model developed in 
this paper can be easily adapted for studying the influence of hot and dense medium which appears 
in ultrarelativistic heavy ion collisions and what could be fairly crucial step in accomplishing the 
objective. 
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